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The Periodogram and its Optical Analogy . 
By Arthur Schuster, F.B.S. 


(Beceived November 23,—Bead December 7, 1905.) 


I have recently applied the periodogram method to the investigation of 
several fluctuating quantities, and the experience thus gained has led me to 
modify slightly the original definition.* Having always laid stress on the 
fact that the periodogram supplies by calculation the transformation which 
the spectroscope instrumentally impresses on a luminous disturbance, I may 
now enter a little more closely into this optical analogy, and thus lead up to 
what I hope will be the final definition. 

Consider a parallel beam of light falling on a grating, the reflected light 
being collected at the focus of an observing telescope in the usual way. For 
simplicity of calculation I assume that the grating considered is of a 
particular type, which, in a former paper, I have called a simple grating. 
Such a grating only gives two spectra of the first order. 

If <j> (Yt -f x) be the velocity at any point of the incident beam, the 
displacement at the focus of the observing telescope isf 

hi cos/3 -p 
2 tt/NVX ’ 


where 


_ r+*N* 

B = I cos nx <b(Yt + x) dx. 
J-JNA 


a) 


In these equations f denotes the focal length of the telescope, h is the 
length of the lines ruled on the grating, l the width of ruled space measured 
at right angles to the lines, N gives the number of lines, and ft the angle 
between the direction of the optic axis of the observing telescope and the 
normal to the grating. For the sake of shortness n is written for 2ir]\. The 
quantity denoted by X is the wave-length of homogeneous light which would 
have its first principal maximum at the focus of the telescope. It may be 
said to be the wave-length towards which the telescope points, and its strict 
definition is given by the relation 

l (sin oc —sin /3) = NX, 
where a is the angle of incidence. 

In order not to complicate needlessly the calculations, I shall assume that 
the resolving power is sufficient to ensure that at any point of the spectrum 


* 4 Cambridge Phil. Soc. Trans./ vol. 18, p. 107. 
f ‘Phil. Mag./ vol. 37, p. 545 (1894). 
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the vibrations are nearly homogeneous; this involves that the average 
squares of the velocities are sensibly equal to the average squares of the 
displacements multiplied by 47 t 2 V 2 /X 2 . The average square of the velocity 
at the focus of the telescope is in that case— 

ll 2 l 2 COS 2 /3p 2 
/ 2 N 2 X 4 

where for K 2 we must put its average value. This expression represents the 
measure of the intensity at the point considered. Its line and surface 
integrals may be called the total linear intensity and the total intensity 
respectively. 

In observing a spectrum, we associate with a particular wave-length all 
the light which lies in a straight line parallel to the rulings of the grating. 
The distribution of light along a vertical line for nearly homogeneous light 
takes place according to the law aC 2 sin 2 a, where a = 7 rhylf\,y being the 
vertical distance. Multiplying by dy, and integrating from minus to plus 
infinity, the total intensity in a vertical line is found to be X/fh when the 
intensity at the central maximum is unity. With the value for the central 
intensity previously found, we now obtain the total linear intensity associated 
with X to be 

hi 2 COS 2 /3~p2 /0\ 

fWX* * ( } 

Changing the variable, the expression for E takes the form 

f # 0 +NA 

cos n(Vt—x)<h {os) doc, 

Xq 

where x 0 is put for Yt— JjSTX. 

Write 

f ^ 0 +E'A pc 0 +NA 

cos nx <f> {x) dx ; B — sin rix <£ (x) dx. 

Xq J Xq 

The mean value of B 2 is then equal to the mean value of 

i(A 2 + B 2 ). 

A grating such as that to which the above equations apply forms two 
spectra and absorbs part of the light; we must now estimate what 
fraction of the incident beam is utilised to form the spectrum under 
consideration. Eor this purpose we imagine homogeneous light to fall on the 
grating, and put cf> (x) = cos nx. The mean value of | (A 2 -f B 2 ) is then 
easily found to be N 2 X 2 . By substitution into (2) we find that the total 
linear intensity in the central line is now 

7^ 2 cos 2 /3-p 2 
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To either side of the principal maximum the intensity varies according to 
the law #~ 2 sin 2 a, where a is equal to cos/9//\, £ representing a 
distance measured at right angles to the spectroscopic line. The total 
energy measured in the focal plane of the telescope is obtained by multi* 
plying (3) with a"* 2 sin 2 otd £, and integrating. This gives £ hi cos (3. 

If the incident light is normal to the grating, its total energy is-JAZ, the 
factor J representing the fact that we have taken the average square of the 
velocity which is half the square of the greatest velocity as the measure of 
intensity. We conclude that | cos /3 is the fraction of light utilised to form 
the spectrum. 

Taking account in (2) of this, we find that the type of spectroscope 
considered estimates the intensity of light passing through its central 
meridian as being 


where for A 2 and B 2 their mean values are to be substituted. 

To obtain the total light within a small angular distance d/3, we must 
multiply by Frf/8; as NcZ\ = l cos /3d/3, we find that the total energy within 
a range d\ is 


£(A 2 +b 2 )& 


If the total energy of the light incident on the grating is unity, the energy 
assigned by the grating to a range dn is therefore finally— 


In the application of the periodogram it is more convenient to take the 
time as the independent variable. Defining therefore— 



A 2 4-1 $ 2 

(4) becomes equal to ——— dfc. 


7 ^ ttNT * 

Leaving out the constant factor, I now define 

S = (A 2 +B 2 )/NT 


to be the ordinate of the periodogram. The definition differs from the 
previous one by the factor NT, which occurs in the denominator instead of its 
square. 

The present definition is not only justified by the close optical analogy 
which has now been formally proved, but also by the resulting convenience. 
I have previously shown that, in the absence of any homogeneous 
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periodicities, the average of A 2 4-B 2 increases in proportion to the time 
interval, NT, which occurs in the limits of the integrals for A and B. It 
follows that for such variations, the ordinate of the periodogram as at 
present defined is independent of the time limits chosen. This is an 
advantage. On the other hand, the former definition gave directly the 
amplitude of the periodic variation when it was of an absolutely homo¬ 
geneous character. For such homogeneous variation the present ordinate 
increases proportionally to the time interval chosen. 

The optical analogy explains the reason of this, and gives its justification. 
When homogeneous light falls on an instrument of definite resolving power 
the light in the central meridian does not by itself give sufficient indication 
of the intensity of the incident light. It is only when correction has been 
made for the lateral spreading that the true intensity can be deduced, the 
correction depending on the resolving power. It is otherwise when the 
spectrum is continuous, for in that case the light lost by lateral spreading 
is replaced by that which properly belongs to the neighbouring wave-lengths. 
Hence, in this case, the intensity in the central meridian is a true measure 
of the intensity of the incident light. 

It need hardly be pointed out how constant use is made of the fact that 
increased resolving power {not increased dispersion) brings out the homo¬ 
geneous lines of a spectrum by increasing their intensity beyond that of the 
continuous background. It is correspondingly one of the principal advantages 
of the periodogram method that it gives a measure of the resolving power 
necessary to isolate a true homogeneous period from the irregular fluctuations. 

Light is thrown on parts of the previous investigation by a formula 
given by Lord Kayleigh for the intensity to be assigned to the homogeneous 
components of a disturbance. If $(%) be the velocity at any point of a linear 
disturbance, so that the total intensity is 

f (c f>(%)} 2 dx , 

J «*00 

Lord Rayleigh shows that the energy to he assigned to a range dn, where 
n = 27r/X, is 

(A 2 +B 2 )cfc/7T, 

in which 

f +oo f+ 00 

cos kv <j> {v) dv, B == I sin tcv $ (v) dv . 

— oo j —oo 

The average intensity spread over a certain length L may be estimated 
by taking and v 0 + L as lower and upper limits of the integrals, and 
averaging the values obtained by a change of v 0 . The energy per unit 
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length would then he found on dividing the expression in (5) by L. We 
arrive in this manner at equation (4). 

I might have confined myself to this simple deduction had I not wished 
to lay stress on the equations for instrumental resolution. This seemed 
all the more desirable because for absolutely homogeneous radiation a 
definition of the periodogram based on the average intensity per unit length 
would fail. If a simple periodicity exists, its amplitude may easily be 
derived from the ordinate of the periodogram, for, if S be that ordinate, 
the amplitude is 2 (S/NT)£. 

In practical applications the function <j>{t) will generally be given for 
successive intervals (a) of the time. The integrals occurring in (5) are then 
replaced by summations, unless a harmonic analyser is used. It is most 
convenient to write in this case 


s = (n—\) a 


/2tt \ 

8 = (n—1) a 

/27T \ 

(— s ) > 

B = 2 sin 

(— s ) 


5 = 0 

\ n ) 


( 6 ) 


S =(A 2 + B 2 )a 2 /NT, 

where <£ s represents the values which <f> (t) takes at the successive times 
considered. 

If we take p to be equal to the total number of separate values of </> ( t) 
used in the calculations, we may put 

S = (A 2 +B 2 )*/p. (7) 

If a harmonic analyser be used, and a b are the two Fourier coefficients, 

S = i (a 2 + b 2 ) NT = i (a 2 +b 2 ) *p, 


where ap represents the complete time interval to which the Fourier analysis 
has been applied. 









